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ABSTRACT 

We present a detailed analysis on the phases of the WMAP foregrounds (synchrotron, free-free and 
dust emission) of the WMAP K-W bands in order to estimate the significance of the variation of the 
spectral indices at different components. We first extract the spectral-index varying signals by assuming 
that the invariant part among different frequency bands have 100% cross-correlation of phases. We then 
use the minimization of variance, which is normally used for extracting the CMB signals, to extract 
the frequency independent signals. Such a common signal in each foreground component could play a 
significant role for any kind of component separation methods, because the methods cannot discriminate 
frequency independent foregrounds and CMB. 
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1. INTRODUCTION 

The release of the first year data from the Wilkinson 
Microwave Anisotropy Probe (WMAP ) provides a unique 
opportunity to investigate the properties of whole-sky cos- 
mic microwave background (CMB) anisotropies at un- 
precedented resolution and sensitivity to date (Bennett et 
al. 2003b; Hinshaw et al. 2003a,b). The WMAP science 
team also produce several whole sky foreground maps 
of synchrotron, free-free and dust emission from all the 
frequency bands at the range 23-94 GHz with unprece- 
dented precision (Bennett et al. 2003c). The analysis of 
the foreground properties is important as it is an indi- 
rect testing of the derived CMB signal properties (the an- 
gular power spectrum and its corresponding optimal val- 
ues of the cosmological parameters for the best fit ACDM 
model (Tegmark, de Oliveira-Costa & Hamilton 2003; de 
Oliveira-Costa, Tegmark, Zaldarriaga & Hamilton 2003; 
Efstathiou 2003), and it also sheds light on the statistical 
properties of the CMB signal, i.e. Gaussianity of the CMB 
anisotropies (Komatsu et al. 2003; Chiang et al. 2003; Park 
2004; Naselsky, Doroshkevich & Verkhodanov 2003, 2004; 
Eriksen et al. 2004a,b; Hansen et al. 2004). 

In the framework of inflation paradigm, the angular 
power spectrum Ci provides all the information about 
the CMB fluctuations when the anisotropies constitute a 
Gaussian random field (Bardeen et al. 1986; Bond & Ef- 
stathiou 1987). Therefore, non-Gaussianity, either with 
primordial origin, or residuals from systematics or fore- 
ground component separation, shall inflict the accuracy of 
the cosmological parameters. Recently there are reports 
on detection of non-Gaussianity from 1-year derived CMB 
maps, which is attributed to foreground residuals or the 
systematic effect correction (Chiang et al. 2003; Naselsky, 
Doroshkevich & Verkhodanov 2003, 2004; Eriksen et al. 
2004b; Hansen et al. 2004). This issue therefore needs to 
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be dealt with care before we reach any conclusions on, 
among other important issues, the statistical properties of 
the CMB anisotropies. 

The WMAP science team use Maximum Entropy 
Method to produce five maps at K-W bands for the syn- 
chrotron component, five maps for the free-free emission. 
They also incorporate extrapolation of the galactic dust 
emission from 100 microns (Finkbeiner, Davis and Schlegel 
1999) to WMAP frequency range. In this paper we exam- 
ine the properties of the WMAP derived foreground maps, 
using an approach based on the phase of spherical har- 
monic coefficients. From theoretical point of view, the 
phases of the major foreground components listed above 
are extremely important. They provide an opportunity to 
develop and to implement so called "blind" method of the 
CMB and foregrounds separation not for the WMAP data 
only, but for the upcoming high-resolution PLANCK data 
as well. 

The CMB and the combined foreground signals: syn- 
chrotron, free-free, dust emission and so on come from the 
same data sets. For the method by the WMAP science 
team for the CMB signal extraction, that by Tegmark, 
de Oliveira-Costa & Hamilton (2003) (TOH hereafter), 
Naselsky et al. (2003), and by Eriksen et al. (2004b), pro- 
nounced correlations of phases (hence correlations in mor- 
phology) between different frequencies at each foreground 
component arc necessary (Tegmark & Efstathiou 1996; 
Bouchet & Gispert 1999; Hobson et al. 1998; Stolyarov et 
al. 2002; Delabrouille, Patanchon & Audit 2003; Tegmark, 
de Oliveira-Costa & Hamilton 2003; Patanchon et al. 2003; 
Barreiro et al. 2003; Martinez- Gonzalez et al. 2003), i.e. 
the foregrounds at different frequencies should only differ 
by amplitudes not by morphology. Using the derived maps 
for the synchrotron, free-free and dust emission from the 
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WMAP web site 5 we demonstrate such correlations (see 
Section 1). Not surprisingly, these phases are not strongly 
correlated, especially for the synchrotron emission. Such 
deviation in morphology is mentioned by WMAP science 
team as the so-called spectral index variation effect (Ben- 
nett et al. 2003c). The aim of this paper is to find the 
foreground components of signal from K-W band which 
cause deviations in morphology in the foregrounds (Sec- 
tion 2 and 3) and some components which do not have 
any frequency dependence (Section 4). The latter creates 
some problems for the standard tools for CMB component 
separation and determines the corresponding error of the 
tools. Below we will mainly discuss the synchrotron, free- 
free and dust components, but our method can be easily 
expanded on all kind of foregrounds, which are important 
for the component separation tools for the Planck mission. 

This paper is arranged as follows: in Section 2 we 
use phase analysis to develop a method for the spectral- 
index varying signals in the foreground maps. For the 
WMAP foreground maps the derived spectral-index vary- 
ing signals are shown in Section 3. And in Section 
4 we apply the minimization of the variance on the 
WMAP foreground maps to extract the frequency inde- 
pendent signals for each foreground component. We use 
circular statistics to test the cross correlations of those 
maps in Section 5. Conclusions and discussions are in 
Section 6. 

2. SEPARATION OF THE SPECTRAL-INDEX VARYING 
SIGNAL IN THE WMAP FOREGROUND MAPS 

We formulate the question as follows: how can we ex- 
tract the common phases for each kind of the foregrounds 
from the WMAP derived foreground maps? 

For that purpose we will apply circular phase statis- 
tics (Fisher 1993), using Lagrangian weighting coefficients 
Wj (£) , (J2j w j {£) = I) ! similar to the TE96 optimization 
(Tegmark & Efstathiou 1996; Tegmark, de Oliveira-Costa 
& Hamilton 2003), but for the phases at each band j. For 
statistical characterization of the temperature fluctuations 
on a sphere we express each foreground signal as a sum 
over spherical harmonics 
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where \ag m \ and 4>e m are the moduli and phases of the 
coefficients of the expansion. 

In practice, we use the Healpix package (Gorski, 
Hivon & Wandelt 1999) to decompose each of the 
WMAP foreground maps at K-W bands for the spheri- 
cal harmonic coefficients ag m . Then we take the whitened 
foreground maps for each frequency band j: 
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and using the weighting coefficients Wj(£) on the phases 
we have the map, for example for the synchrotron emis- 
sion S£ m ■ 



exp (i$ 



= n^ (Mir)^Uxp ^E^W^' 5 j • (3) 
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The weighting coefficients are designed to produce the map 
whose phases $f m = J2j w j(^)^hn' S do not depend on fre- 
quency. The map Se m represents the common part of the 
foregrounds (in morphology) among the K-W bands. Note 
that we can reach different sets of weighting coefficients for 
the maps of the free- free {Ttm) an d the dust (T>i m ) emis- 
sions. 

2.1. Determination of the weighting coefficients 

To obtain the weighting coefficients Wj(£) it is natural 
to assume that 3>f m is orthogonal to the phase differences 
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1 — 5 corresponds to the K, 



Ka, Q, V and W bands, respectively. Taking into account 
of J2j w j(?-) = I in Eq.(3) we have (for generality we drop 
the superscript <S) 



where = - The orthogonality of and 
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However, considering 2tt periodicity of the phases, a sim- 
ple generalization of Eq.(5) should be considered, which 
can be obtained by the following way (Fisher 1993). We 
introduce the trigonometric variables 



x lm = sin$ £m ; Vtil = sillD tm (6) 
and the cross-correlation coefficients r defined as follows 
should produce zero value due to orthogonality, i.e. 
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r{£) = Esin^sin^ = 0. 
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As one can see from Eq.(7), if the phase difference and 
&i m are small (<C 7r/2), we will have the same equation 

for the weighting coefficients as Eq.(5). If is small, 
but <&e m is comparable with tt/2 we will have non- linear 
equation for the weighting coefficients that can be solved 
numerically 

2.2. Spectral-index varying foregrounds 

We now can derive, from the given coefficients w, the 
maps for each component of the foregrounds (synchrotron, 
free-free and dust emission) with the common phases of 
components among all frequency bands as follows: 
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where j denotes the frequency bands, S, T, T> denote the 
synchrotron, free-free and dust emission, respectively, and 
®tm corresponds to each component at each band. The 
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maps of spectral- index varying signal at each band for each 
foreground component are 



L lrn 



exp 



(9) 




Taking into account Eq.(9) one gets the map of fluctuation 



tra 



where = 



(10) 

$£ m . The map from Eq.(10) repre- 
sents all variation of the phase for all m > 1 harmonics of 
the signal, if A^fj <C 1. As one can see, for that map the 



moduli are la 



I'm 



and the phases are the initial 
phases rotated by the angle tt/2, if > and by the 
angle 3tt/2, if Ag\ < 0. 






3. THE SPECTRAL-INDEX VARYING SIGNALS IN THE 
WMAP SYNCHROTRON, FREE-FREE AND DUST MAPS 

We implement the method we have developed to extract 
the peculiar signals for the synchrotron, free-free and dust 
emission for the WMAP frequency bands. We do not in- 
clude the K band in our analysis due to the same reason 
as the WMAP science team mentions. The maps for the 
spectral-index varying signals are grouped in three: syn- 
chrotron, free-free and dust emission. 

In Fig.l we plot the maps for the spectral- index varying 
signal of the synchrotron emission for the WMAP Ka-W 
bands. The m = modes in these spectral- index vary- 
ing maps have zero amplitude. It is the reason why syn- 
chrotron maps looks peculiar. 

Fig. 2 and 3 are the spectral-index varying signals for the 
free-free and dust emission, respectively. Again the m = 
modes are naturallly excluded. 



Fig. 1. — The derived spectral-index varying signals for the syn- 
chrotron component. From top to bottom are Ka, Q, V and W 
bands. 

For the WMAP dust maps the spectral-index varying 
signals are relatively weak in comparison with those from 
the synchrotron and free-free maps. However, it is in or- 
der of 30% of the amplitude of ILC signal which should be 
taken into account seriously. 

In Fig. 4 we display the spectral-index varying signals, 
the sum of the synchrotron, free-free and dust emission. 




Fig. 2. — The derived spectral-index varying signals for the free- 
free emission. From top to bottom are Ka, Q, V and W bands. 



Fig. 3. — The derived spectral-index varying signals for the dust 
emission. From top to bottom are Ka, Q, V and W bands. 
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Fig. 4. — The derived spectral-index varying signals of the sum of 
the synchrotron, free-free and dust emissions. From top to bottom 
are ka, Q, V and W bands. 



4. THE FREQUENCY INDEPENDENT SIGNALS (FIS) IN 
THE WMAP FOREGROUND MAPS 

In this Section we use the minimizatoin of the variace 
in a similar fashion on the WMAP foreground maps for 
the frequency independent signals (hereafter FIS). This 
method is mainly used by the WMAP science team to pro- 
duce the CMB map (Bennett et al. 2003c) using a weighted 
combination of the 5 bands outside the Galactic plane. 



Tegmark, de Oliveira-Costa & Hamilton (2003) (hereafter 
TOH) perform an independent foreground analysis from 
the WMAP data using weighting coefficients dependent 
not only on the angular scales (£) but also on the Galactic 
latitudes. Below we use the TOH method for the syn- 
chrotron, free-free and dust maps in order to extract the 
frequency independent component (or partly independent 
for some of the frequency range) without any assumptions 
about the statistical nature of the derived signals. More- 
over, we use the whole sky for analyses and do not apply 
any Galaxy plane cut-off and masks, nor do we dissect the 
whole sky into disjoint regions. 

Using minimization scheme by TOH and Tegmark & Ef- 
stathiou (1996), we implement the blind method for the 
FIS for synchrotron, free-free and dust emission which 
have the same phases for K-W bands. Unlike Section 
1, here we are looking for some common component of 
the foregrounds, which is frequency independent at the 
WMAP frequency range. Moreover, such a common com- 
ponent could play a significant role for any kind of compo- 
nent separation methods, because the methods cannot dis- 
criminate frequency independent foregrounds and CMB. 

4.1. The synchrotron FIS maps 

We firstly group the maps in the following combinations 
for the WMAP synchrotron maps: Ka-Q, Ka-V and Q- 
V. For each pair of the maps (Ka-Q, Ka-V and Q-V) 
we want to reach the common Synchrotron emission map 
(af m = S) which has the same phases for each pair of the 
bands and for each £, m mode. We use the optimization 

coefficients T {l) = , similar to that by TE96 Tegmark 
& Efstathiou (1996), TOH, Naselsky et al. (2003): 



2 2 

-.(<)„(<) =V , r(i)Ji) 



S = af m = £rWa«^£rWa(' 



(11) 



i=l 



where i = 1, 2 correspond to the two bands in a pair. The 
coefficients are subject to the following constraints 
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and 



£r« = i, 



= o, 



(12) 
(13) 



where S/ST^ denotes functional derivatives. 
The rW then have the following forms 

(1) _ EL = ^{|a (2) l[l a(2) |-| a(1) l^ (2) ^ (1) )]} 
£ f m= _ f |a (1) -a (2) | 2 ' 

r( 2) _ SU-^la^HlaWl-laWlcoB^W-^ 1 ))]} 



-a< 2 )| 2 



(14) 



where are the phases of the synchrotron signal aW at 
the i-th map. 

Then we obtain the common signal S and the frequency 
dependent part at each j-th bands. 



F£ ) =a( 2 )-S = -r( 1 )(aW-a( 2 )). 



(15) 
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Fig. 5. — The power spectra of the common component from 
the synchrotron emission. Top thick dash line corresponds to the 
WMAP best fit ACDM model. Boxes reproduce the ILC power 
spectrum. Thin long dash line is the common component from Ka- 
Q synchrotron maps. Solid line is the common component from 
Ka— V synchrotron maps. Dash-dot line is that from Q-V and Dash 
dot-dot line is that from V-W. 

In Fig. 5 we plot the power spectra for the FIS taken from 
different combinations of the WMAP synchcrotron maps. 
As one can see, for Q-V bands and Ka-V bands we have 
the frequency independent component. 




Fig. 6. — The FIS from Q-V synchrotron maps. 

In Fig. 6 we show the FIS from Q-V synchrotron maps, 
which mainly concentrated in the Galactic plane, but some 
of the point source residues seen above and below it. The 
power spectra of such FIS shown in Fig. 5 is smaller than 
the ILC power, as shown in Fig. 5. However, due to the 
source clustering in the Galactic plane it can significantly 
perturb the CMB map at Galactic plane region (see for 
example, Naselsky et al. (2003)). 
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Fig. 7. — Reconstructed phases of the FIS from synchrotron 
maps with different combination of the bands. From top to bottom 
panels are cross-correlation between combination of Ka-V (y-axis) 
against Q— V (x-axis), Ka-Q against Q— V, V-W against Q-V, and 
ILC against Q— V. 
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In Fig. 7 we plot the cross-correlations of phases for the 
synchrotron FIS obtained at different frequency bands. 
This Figure confirms that the FIS is common for Ka-V 
bands. 




fju.': Y/ 



4.2. The dust FIS maps 

We show the dust FIS, using the same separation tech- 
nique as for synchrotron emission. 
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Fig. 8. — The power spectra of the common component from the 
dust emission. Top thick dash line corresponds to the WMAP best 
fit ACDM model. Boxes reproduce the ILC power spectrum. Thin 
long dash line is the common component from Ka-Q synchrotron 
maps. Solid line is the common component from Ka-V synchrotron 
maps. Dash-dot line is that from Q-V and Dash dot-dot line is that 
from V-W. 
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Fig. 9.— The FIS from Q-V dust maps. 

Fig. 8 is the power spectra. As shown in Fig. 9 the FIS is 
mainly related to the point sources. Contamination of the 
Galactic plane signal is week. In Fig. 10 we plot the cross 
correlation of phases for the FIS at different combinations 
of the bands. These correlations are very strong for all 
Ka-V bands. 
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Fig. 10. — Reconstructed phases of the FIS from the dust maps 
with different combination of the bands. From top to bottom panels 
are cross-correlation between combination of Ka— V (y-axis) against 
Q— V (ir-axis), Ka-Q against Q-V, V-W against Q— V, and ILC 
against Q— V. 
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Fig. 11. — Power spectra for the common component from free- 
free emission. The definition of the lines are the same as in Fig. 5 

The power spectra arc different for different frequency 
bands and phases of the signal are not highly cor- 
related. Actually, the free-free component of the 
WMAP foreground is the only one which does not have 
FIS. 

4.3. The free-free FIS maps 

For the free-free emission the derived FIS is extremely 
weak, as one can see from the Fig.ll. 

5. CROSS CORRELATION OF PHASES BETWEEN THE 
DERIVED CMB MAPS AND THE FOREGROUNDS 

The aim of this Section is to show how the variation of 
the spectral index in the foreground components can man- 
ifest themselves in the CMB signal, which should be clean 
from the foreground contamination. We use the whole-sky 
maps derived from the WMAP data using different sepa- 
ration methods. Namely we use the WMAP Internal Lin- 
ear Combination (ILC) map, the foreground cleaned map 
(FCM) and Wiener filtered map (WFM) by TOH, and the 
phase cleaned map (PCM) by Naselsky et al. (2003). 

We use the trigonometric moment statistics to counter 
the circular nature of phases (Fisher (1993), see also Nasel- 
sky, Doroshkevich & Vcrkhodanov (2003)). We define the 
following trigonometric moments: 



C(t,At) 
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m— 1 



t-t 



r 2 {t, At) = [C 2 (t, At) + S 2 (t, At)] t, 



(16) 



where At — £ + t and ^e m are the phases of the 
foreground-cleaned signals (ILC, FCM, WFM, PCM) at 
each tm harmonic and <&e m are the phases of the fore- 
grounds. As we have shown, the most dangerous compo- 
nent of the foreground at each frequency band is the spec- 
tral index varying tail in the synchrotron, free-free and 
dust emission, shown in Fig. Below we plot the cross cor- 
relation function S, C and r 2 for each of the foreground 
cleaned signals. Using the statistics we look for signatures 
that deviate from the 68% CL. If the phases of the cleaned 



map arc highly correlated with those of the foregrounds, 
the C statistic has asymptotics to unity, whereas the S 
to zero. For the case of non-correlations, both C and S 
statistics have t 1 ^ 2 asymptotics due to "random walk" . 

In Fig. 12 and 13 we plot the S, C and r 2 statistics on 
the foreground-cleaned and the foreground maps. The top 
3 rows of Fig. 12 is the statistics for ILC and foregrounds, 
in which from left to right columns are for At — 0, 1 
and 2, respectively. One can see significant cross corre- 
lations at multipoles I = 6, 8, 14-16, 21, 39 and 47 of 
ILC with Q, V and W foregrounds at 95% CL and at 
t = 11, 17, 26, 36 and 47 at 99% CL. The bottom 3 rows 
of Fig. 12 is the statistics on PCM and the foregrounds. 
PCM is produced without galactic cut-off and any prelim- 
inary information about foreground properties (Naselsky 
et al. 2003). The possible non-Gaussianity in the PCM is 
the residues of bright point sources in the Galactic plane 
region covered by WMAP Kp2 mask. These residues have 
pronounced cross correlations with the synchrotron com- 
ponent as shown in Fig. 6. In Fig. 13 on the top 3 rows we 
can also see clearly some multipoles with cross correlations 
above 2— a and for t = 30 the correlation is above 99% CL. 

To show the possible multipole range in which con- 
tamination of non-Gaussian features can be important, 
in Fig. 14 wc plot the difference between the best-fit 
WMAP ACDM, ILC and TOH FCM power spectra, nor- 
malized to the ACDM power. The intriguing features ap- 
pear at multipoles I = 2, 5, 7, 20 — 23, 40, in which the 
peaks of ACe/Ce lies above 68% CL. 
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Fig. 14.— The difference between the best-fit WMAP ACDM, 
ILC and TOH FCM power spectra, normalized to the ACDM power. 
The diamond shape curves represent the cosmic variance. 



6. CONCLUSIONS AND DISCUSSIONS 

In this paper we raise the issue of how different the 
foreground maps at different frequency bands. From each 
WMAP foreground component: synchrotron, free-free and 
dust maps, we look for spectral-index varying signals. By 
assuming that the spectral-index invariant part has 100% 
phase correlation (hence the same morphology) among the 
different frequency bands, we develop a new method for 
extracting the spectral-index varying signals, which corre- 
spond to perturbations in phases. We have shown that this 
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Fig. 12. — The circular statistics for the cross correlation of phases between the ILC and the foreground maps (top 3 rows) and between 
the PCM and the foreground maps (lower 3 rows). Each 3 rows (downwards) are S, C and r 2 statistics, respectively. The 3 columns from 
left to right are for A£ = 0, 1 and 2, respectively. In each panel, the star sign marks ILC-Ka foreground cross correlation, the plus sign marks 
ILC-Q, the diamond sign marks ILC-V and the box sign ILC-W. The solid curves denote I— a CL and the shaded areas cover 1000 Gaussian 
realizations from Monte Carlo simulations. The right column are S, C and r 2 statistics between the PCM and foregrounds for A£ = 0, 1 and 
2, respectively. 
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Fig. 13. — The circular statistics for the cross correlation of phases between the TOH FCM and foreground maps (top 3 rows) and between 
the WFM and foreground maps (lower 3 rows). The notation is the same as Fig. 12. 
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varying component has a characteristic in phases: turning 
7r/2 or 37r/2 from the highly correlated foregrounds. 

We further use the standard minimization of the vari- 
ance to extract the frequency independent signals from all 
the frequency bands. This component can play a crucial 
role for the current blind method for component separa- 
tion in CMB studies. Our method is not only useful for 
the foreground analysis on the WMAP data, it shall be 
very crucial for the upcoming Planck data analysis. 
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